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FuU one-loop radiative corrections are calculated for low energy neutrino- 
nucleon quasi-elastic scattering, Pe+P — *■ e^-\-n, which involves both Fermi 
and Gamow- Teller transitions, in the static limit of nucleons. The calcula- 
tion is performed for both angular independent and dependent parts. We 
separate the corrections into the outer and inner parts a la Sirlin. The 
outer part is infrared and ultraviolet finite, and involves the positron ve- 
locity. The calculation of the outer part is straightforward, but that of 
the inner part requires a scrutiny concerning the continuation of the long- 
distance hadronic calculation to the short-distance quark treatment and 
the dependence on the model of hadron structure. We show that the log- 
arithmically divergent parts do not depend on the structure of hadrons 
not only for the Fermi part, but also for the Gamow- Teller part. This 
observation enables us to deal with the inner part for the Gamow- Teller 
transition nearly parallel to that for the Fermi transition. The inner part 
is universal to weak charged-current processes and can be absorbed into 
the modification of the coupling constants up to the order of the inverse 
proton mass 0{l/mp). The resulting 0{a) corrections to the differential 
cross section take the form [1 + 6out{E)] [(1 + C)(l>^ + + , 
where (1) and (cr) stand for the Fermi and Gamow- Teller matrix elements; 
the outer correction 6out is positron energy (E) dependent and takes dif- 
ferent functional forms for the angular independent and dependent parts. 
All factors are explicitly evaluated. 



(1) 



2 



TEXT PRINTED ON FEBRUARY 1, 2008 



PACS numbers: 12.15.Lk, 13.40.Ks, 13.15.+g 

1. Introduction 

The purpose of this paper is to present a calculation of the radiative 
correction to neutrino-nucleon quasielastic scattering (inverse beta decay), 

i>e+P — >-e+ + n, (1) 

whose experimental accuracy now approaches the level that makes radiative 
corrections a practically important problem. The calculation is similar to 
that for neutron beta decay, which has a long history of no less than forty 
years [1]-[11]. 

The calculation of radiative corrections to neutron beta decay is not 
quite straightforward. There are a number of subtleties. They originate 
from the fact that, at low energies, one has to deal with the proton and 
neutron in four- Fermi theory while the loop integral is divergent: this di- 
vergence is made finite in the electroweak theory, but here we must treat 
quarks and continue the quark calculation to hadronic theory at an inter- 
mediate energy scale. 

To clarify the issues and the scope of the present paper, we start with a 
brief survey of the work done in the past. The pioneering work of Kinoshita 
and Sirlin [1], long before the development of renormalizable gauge theories, 
evaluated photon exchange corrections to the neutron beta-decay within 
four-Fermi V—A theory of weak interactions. They showed the cancellation 
of infrared divergences and derived the correct velocity dependence of the 
final electron, taming the ultraviolet divergence with cut-off theory. 

Repeating the calculation of Kinoshita and Sirlin, Berman and Sirlin 
[2] speculated that the logarithmic divergences are not affected by strong 
interactions. This was later proven by Abcrs ct al. [3]. The theorem that 
the logarithmic divergences are universal is based on the conserved vector 
current with the use of the current algebra technique [12] [13]. This theorem, 
however, applies only to the logarithmic divergent part of super-allowed 
Fermi transitions. More precisely speaking, it applies only to the purely 
vector-current contribution in the Fermi transitions: it does not apply to 
the logarithmic divergent part arising from the axial current which appears 
on the loop level by interference. The divergences due to the axial current 
(which are linear in the axial- vector coupling g^) in general depend on the 
model of hadron structure and strong interaction. Abers et al. [3] proposed 
a model (Ai-exchangc model) to evaluate such contributions. 

Sirlin [4] has gone one step further to separate radiative corrections of 
neutron beta decay into outer and inner parts. The outer part corrections 
are both infrared and ultraviolet finite, gauge-independent, and contain full 
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electron's velocity dependence, thus depending on specific nuclei that un- 
dergo beta decay. The outer part is not affected by strong interactions. The 
inner part is independent of electron's velocity, infrared finite but contains 
the ultraviolet divergences. This separation enables one to absorb the inner 
part into a universal multiplicative correction factor of the vector coupling 
constant that is relevant to all beta decays. 

Electroweak gauge theory renders the ultraviolet divergences of the inner 
part finite. The calculation [7, 8] conventionally divides the integration 
region of the virtual gauge bosons into long- and short-distance parts: 

(i) < \k\^ < M^, (ii) < |/c|2 < GO, (2) 

where k is the (Wick-rotated Euclidean) momentum of the virtual gauge 

bosons, and the mass scale M, introduced by hand, divides the low- and 
high-energy regimes and is supposed to lie between the proton-neutron 
masses {nip and m^) and the {W,Z) boson masses, mw and mz- Old- 
fashioned four-Fcrmi interactions are applied to the proton and neutron in 
region (i) , and the mass scale M is regarded as the ultraviolet cutoff of the 
QED (i.e., purely photonic) correction. In region (ii), electroweak theory is 
used for quarks and leptons, and M is a mass scale that describes the onset 
of the asymptotic behaviour. The concern is whether the results in (i) and 
(ii) join together smoothly. In fact, the theorem mentioned above guaran- 
tees that the logarithmic terms coming from the pure vector current in the 
super-allowed Fermi transitions have the common coefficient in (i) and (ii), 
and the integrals continue smoothly; thereby M-dependence disappears^. 

Radiative corrections proportional to fvQA, on the other hand, requires 
some intricate treatment. (Here /y = 1 denotes the vector coupling con- 
stant parallel to gA', we put fy to trace the vector current contributions.) 
One conceivable way proposed by Sirlin [7] to deal with the axial current 
contribution is to use the free quark model for the short-distance contribu- 
tion in the asymptotic regime, and to work explicitly with the proton and 
neutron for the long-distance part, by introducing nucleon's electromagnetic 
and weak form factors to render the logarithmic divergence milder; the mass 
scale M then remains as a lower cut-off for the asymptotic regime in the 
region (ii). This calculation is necessarily model-dependent. Marciano and 
Sirlin [9] (see also Towner [10]) evaluated the axial current contribution this 
way, and their calculation has widely been employed in 0+ — O"*" beta decay 
phenomenology. 

All the development after Kinoshita and Sirlin, however, is restricted to 
Fermi transitions, whose most important application is to 0"*" — *■ 0+ beta 



^ Contrary to the logarithmically divergent part that is universal, the constant part 
may in principle receive correction from strong interactions. 
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decays; little progress has been made regarding Gamow- Teller transitions, 
except for the early studies which showed that the outer correction of Sirlin 
is also applicable to Gamow-Teller transitions. 

The full calculation of the radiative correction to (1) requires not only 
the treatment of the Fermi transition, but also of the Gamow-Teller transi- 
tion. QED radiative corrections to (1) were already calculated by Vogel [14] 
and Fayans [15] within the cut-off theory with the point nucleons. From 
these calculations one can extract the outer part whereas the inner part is 
basically left untouched, thus circumventing the problems concerning the 
axial current complications if the inner part is empirically evaluated from 
neutron life time. These calculations, that lack the evaluation of the inner 
part, however, would not completely elucidate the structure of radiative 
corrections in quasi-elastic neutrino scattering. 

In this paper wc attempt to calculate full one-loop radiative corrections 
to (1) for both angular independent and dependent parts, including those 
to the Gamow-Teller transition. We prove that the coefficient of the log- 
arithmic divergences arising from the purely axial vector current are not 
affected by strong interactions, which is parallel to the theorem for the 
Fermi transition^. The constant term may receive extra contributions from 
non-conservation of the axial current, but our observation opens up a way 
to combine calculations in two regions in (2) smoothly for the Gamow-Teller 
part. 

The vector current contribution to the Gamow-Teller transition is model 
dependent, just as much as the axial current part that contributes to the 
Fermi transition. Our treatments will be reciprocal between the Fermi and 
the Gamow-Teller transitions. 

Our observations allow the full calculation of the radiative correction to 
the Gamow-Teller transition nearly parallel to that to the Fermi transition, 
though leaving somewhat more room for extra contributions from strong 
interactions to the constant term. We present the full expression of the 
outer and inner radiative corrections to neutrino-nucleon scattering in the 
static limit of nucleons. The application to other charged current processes 
is straightforward. 



^ This proof is implicitly included in the work of Sirlin [16], which discussed the univer- 
sality of the logarithmic divergence for the semi-leptonic processes in the electroweak 
theory. We thank Professor Sirlin to attract our attention to his work. 
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2. Preliminaries 

We use the four-Fermi theory to derive radiative corrections in region 
(i). The matrix element of the Born term for the process (1) is given by 



V2 



V,{jPu)l^{l-l^)Ve{t) [Un{jPl)Wx{j>l,P2)Up{jP2)] , (3) 



where Gy = GpCOsOc with the universal Fermi coupling Gp and the Cabibbo 

angle 6c- The spinors of the neutron, proton, positron and antineutrino 
are denoted by Un, Up, Ve, and v,y, respectively, with momenta specified in 
parentheses. 

The vertex of the hadronic weak current in general depends on micleon's 
momenta, thus written as Wx{pi,p2), but is well-approximated at low en- 
ergies by 

Wxipi,P2) = l\{fv - QAl^), (4) 

where /y = 1 and qa = 1.2670. We retain the constant /y to trace the 
vector current contributions. We do not consider the terms of 0{l/mp). 
The differential cross section is given in terms of the invariant amplitude 

M, 

dajpe +P — > e+ + n) ^ 1 E(3 v^| w|2 
^(0086*) 647r mpTrinEi, ^ 

= ^ E-'a {A{f3) + BWcosO} , (5) 

where 9 is the angle between incident antineutrino momentum and 
positron's, The spin summation is taken over all external fermions. The 
energy of the positron in the final state is denoted hj E = Ei, + nip — nin, 
El, being the incident antineutrino energy; (3 = sjE'^ — ml/ E is the velocity 
of the positron. The forward-backward asymmetry is given by 

Upon spin summation for the Born amplitude (3), we arrive at 
^ = 2,2GlmnmpEE, {(/^ + 2,g\) + {fl - g\)f3co^e] . (7) 

spin 

In the tree level, 

A{P)=Ao = f^ + Sgl, B{P) = Bo = f^-g\. (8) 

We will write fy = fvi^y, and 3g\ = g\{cr)'^ to make the Fermi and 
Gamow-Teller contributions explicit. 
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3. QED corrections 

The QED radiative corrections in A{I3) were partly obtained by Vogel 
[14] and Fayans [15]. These authors, however, did not calculate the inner 
part correction. In this paper we compute the full corrections to both A(P) 
and B(P). Some of our formulae presented in this paper are already obtained 
in [14] and [15]. Nonetheless, since the derivation is not given in the above 
references, we give sketches of calculations so that interested readers can 
check them easily. 

The diagrams we consider are depicted in Figure 1, where (v) is the 
vertex correction, (s) is the self energy correction and (b) is bremsstrahlung. 
We consider the static limit for nucleons, = {pi — p2)^ <C rrip. Our 
calculation is done in the Feynman gauge throughout. 



iPi) (^) 




(b) 

Fig. 1. QED corrections: (v) vertex correction, (s) self-energy correction, and (b) 
bremsstrahlung. 
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3.1. Vertex corrections 
The vertex correction is given by 

ct^k 1 



1 



1 



\/2 7 {2ttY {k- ly -ml {p2-kf -ml k'^ 
xv,-f^{l - r') {^■{k-i) + me} -r^veii) 
xun{pi)Wx{pi,P2 -k){'j- {p2 -k)+ mp} -1^Up{p2), 



(9) 



where A is the photon mass to regulate the infrared divergence. 

Using the identities obtained with the use of the Dirac equation for the 
positron and proton, 

{^■{k-t)+me]l>'ve{l) = {{k-2£r+ia^'-k,}ve{£), (10) 
{7 • (P2 - k)+ mp} 7^np(p2) = {(2p2 - k)f, - ia^^k'"} Up{p2), (11) 



we decompose (9) into three parts (see [4]), 



(12) 



Here M^'"^^ picks up the product of {k — 21Y in (10) and (2p2 — k)^ in 

(11), and at the same time W\{pi,p2 — k) is replaced by W\{pi^p2). It has 
apparently the same gamma matrix structure as the Born term (3), and is 
then written as 



m 



e^I{f3) X 



"I 



d^k 



(k - 21) ■ {2p2 - k) 



(27r)4 _ £)2 _ - kf - ml] {k"^ - A^} 

The integral /(/3) in (14) is given in Appendix of [3]; the real part reads 



(13) 
(14) 



'M2 



1 + log ( ^ ) - 9tanh- V log ( ^ 



mi 



A2 



2/3 
1 + /3 



--(tanh-i/3)2 



where 



L{z) = ^log(l - t) 



(15) 



(16) 



is the Spence function. The ultraviolet cutoff M has been brought about 
by the cut-off of the integral in region (i) of (2). 
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The term represents the combination of ia^'^ky in (10) and {2p2 — 

k)fj, in (11), and W\{pi,p2 — k) is again replaced with W\{pi,p2), i.e., 

^^(^2) ^ ±^ 2 f 1 1 

V2 J (27r)4 (k - £)2 - ml {p2 - -mjk'^- 
xv,,^^{l--i^)m^'''k,ve{l) 

XUn{pi)Wx{pi,P2){'ip2 - k)f,Up{p2). (17) 

The gamma matrix structure of the nucleon currents in (17) is the same as 
that in (3), although that of leptons is not. The contributions of M.'^'"^^ and 
are thus evaluated without referring to details of the hadronic part of 
the currents (sec Appendix A). The radiative correction to the cross section 
is given as the tree-one-loop interference amplitude, which takes the form 

^ {(>f("^) +>f(''^))>f(°)* + (Af^"^)* +>f(''^)*)M(°)} 



spin 



= 32GvmnmpEEi, 



im + /(/?)*) + ^/3tanh-i/3| (/^ + Sgl) 



+ je^ (liP) + liP)*) + ^itanh-i/sj (/^ - g\)Pcose . (18) 
The remaining terms are collected in which reads 

>t(-3) ^ ±(. 2 \ \ 

^/2 J (27r)4 {k - if - ml {p2 - kf -mlk^- 
xv,j^{l - 7^) {{k - 2iY + ia^''ky}vS) 

XUn{pi)Rij.x{pi,P2,k)Up{p2), (19) 

where 

Rixx{Pi,P2,k) = {Wx{pi,p2-k)-Wx{pi,P2)}{2p2-k)fj, 

-iWx{pi,P2-k)a^,k'' (20) 
~ -iWx{pi,P2-k)a^,k\ (21) 

in the approximation of the point nucleons. It is only this Al'-""^-' term that 
depends on the details of the structure of the weak currents. It is clear from 
the powers of k that this term is infrared convergent. The explicit use of the 
four-Fermi interaction (3) gives an ultraviolet divergence. A straightforward 
calculation (see Appendix B) shows that 
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spin 



= 32GymnmpEEiy 



(22) 



where 



GT 



e 

8^ 



e 

8^ 



3 
4 



5i I ^log fv9A { ^log 1^1+7 



7 



mi 



5 



(23) 



(24) 



The superscripts F and GT refer to the contribution to the Fermi and 
Gamow- Teller transitions. Note that the correction from AK^-^) is written as 
multiplicative factors on the coupling constants for both Fermi and Gamow- 
Teller parts, without disturbing the original angular dependence. 

3. 2. Self energy corrections 

After subtracting the pole term by mass renormalization, we are left 
with the wave function renormalization. The wave function renormalization 
constant for the fermion of mass m is given by [17] 

for the on-shell renormalization. The self energy correction is given by 

= |y^K)-l + ^/^K)-l|7WW, (26) 
and the contribution to the cross section is 

2{^Z2(me)- 1 + ^Z^imp) - l| l^^^f- (27) 



spin 



spin 
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We see that the logM^ dependence in (18) is cancehed by that in (27), 
i.e., the correction from M^'"^'^ + M.^'"^'^ + Al^*^ is ultraviolet finite. The 
infrared divergence remains that is cancelled after the bremsstrahlung con- 
tribution A^(^) is added. The combination M'^'"^^ + TW^''^) ^ j^(s) j^{b) 
is identified as the outer part radiative correction of Sirlin (up to a numeri- 
cal constant) for the neutrino-nucleon quasi-elastic scattering. It is easy to 
show that this term is gauge independent. On the other hand, still 
contains ultraviolet divergence, and it also receives hadronic complications. 
This contribution is taken to be a part of the inner correction. 



3.3. Bremsstrahlung 
The cross section of single photon emission in Figure 1(b) is given by 



a{ve + P — e"*" -I- n + 7) 



(2 



1 1 r d-'e r 

7r)5 ' SmpTUnE^ J 2^ J 



2u 



2 ^ 

spin 



(28) 



where E' = ^ d"^ -\- and cu = Vk^ + X^. In the nucleon static limit, 
find 



we 



(fe)|2 



spin 



GAmpmnGye^ 



[{k + e'f - ml] 



{f^ + 3gl)E, 



E(i--^-^]Hk-i> 



Hfv-9\) 



{£' + k)-pj£ 



.l2 



i/\2 



+ {k-£'){k-p,) 



{k-t){k-p. 



(29) 



After integrating over photon phase space, the terms proportional to { fy + 
3(7^) in (29) contribute to A{(5), and those proportional to {fy — g\) to 
B{P). 

Putting (29) into (28) and integrating over photon momentum, we obtain 

da{i'e + p — > e+ -I- n -I- 7) 
d{cos9) 

= f^E' (^\p{gm^0)Ao+g^'^\p)BoPcose} , (30) 
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where the function g^'^^\(3), which was calculated by Vogel [14] (partly nu- 
merically) and by Fayans [15] using SCHOONSHIP, is written as 



g(''\p) = 41og 



2{E - me 



A 



1 



tanh"^/3 - 1 



2{E + me 



+4 ( ^ tanh"^ (3-1] log 

VP / V me 

--(tanh-i/3)2 + _^tanh-i/3 + -. 
The other function g^'^'^^ is calculated to give 



2/3 



(31) 



/32 ' 2 IP 
- (4 + 1) ta^-/, + (-^ J + 1) (t<mh- W 



1^ 



2/3 
1 + /3 



1 + /3, 



+4 (^1 - -^tanh- log I 



2me 



1 + 



1\ VT^+VT^' 



(32) 

The treatment of the infrared divergence for g^^^^ is somewhat subtle. Steps 
of the bremsstrahlung calculation are presented in Appendix C. 

3.4- QED Corrections: Summary 
We summarise the QED radiative corrections as 



1 + ^OUt + 



1 



1 + <^out + 



-.{crfgl 



GT' 



m 

GT I xGT'' 



(33) 
(34) 



Here 



"out — "out 



= 5out(^) 

= e2{/(/3) + /(/3)*} + ^/3tanh-i/3 
+2 |^Z2(me) - 1 + ^Z2(mp) - l| 



(35) 
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e 

8^ 




(36) 



is the so-called outer correction, which assembles the contributions from 
M^'"^\ M'''"'^\ M^^\ and M^'^\ This term does not contain infrared or 
ultraviolet divergence, and docs not receive complications due to hadronic 
structure and strong interaction. All positron- velocity dependences are con- 
tained in this correction. We note that the last term 3/4 in (35) is added 
(then subtracted from Sin) to make the definitions of (5out and consistent 
with those of Sirlin for beta decay. The outer correction is common to the 
Fermi and Gamow-Teller transitions, and it can be factored out to the order 
of 0(e2). 

The inner correction (5in' (prime means the correction from QED only) 
arises from Ai ^'"^^ . We find 



s9^' = 




e 
8^ 



3, ^ /vf3, 




(37) 
(38) 



Here the constant 3/4 is subtracted from <1>^ and <I>^^ to define Sf^^' and 6^^' . 
The inner correction is infrared finite, but contains ultraviolet divergences, 
log M^, which are made finite only with electroweak theory. These terms do 
not contain positron- velocity dependent factors, and thus absorbed into the 
vector and axial-vector coupling constants. They also agree with the inner 
correction for beta decay. The inner correction for the Fermi transition has 
been identified in [4]. The first logarithmic divergent parts, both in (5? and 
, do not receive corrections (see Sect. 5). However, the constant terms 
may, in principle, receive corrections from strong interactions. We shall 
argue that effects of the hadron structure are essential for terms proportional 
to qa/ fv in (37) and fv/gA in (38), including the coefficients of logarithmic 
divergences. 

For the correction of the angular dependent part, we have 



"out — "out 



= e2{I(/3) + I(/3)*} + ^itanh-i/3 
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+2 jy^Z^K) - 1 + ^z^K) - 1} + ^ [g^'^KP) + 



e 

8^ 



3 3, m. 



m 



'1-/3 



4_ _ 4\/r^^ 



+4 1 - -tanh" V log - 1 + - 



1 



i\ yrT73 + ^/r^ 



/3 



+ ( - - 4 ) tanh-V + ( ^ - o - ) (tanh-^/?) 



/3 



2/32 



(39) 



The inner correction 5f^ and (^J^'^' are, as evident from (22), identical with 
those that appear in the angular independent correction (33); hence we have 
not attached tildes. 



4. Electroweak corrections and the continuation to low energy 

The short distance correction from the integration region (ii) in (2) is 
evaluated using electroweak theory [7]. For the application of electroweak 
theory, we must deal with quarks, so that the continuation of quark theory 
to low-energy effective theory for nucleons is essential. When wc consider 
corrections relative to muon decay, we only need to consider the photon and 
Z exchange diagrams shown in Figure 2. 




Fig. 2. Photon and Z-boson exchange diagrai 



For a short distance calculation, we can neglect mass and momenta of 
the external fermions. This calculation was done in [7]. A calculation for 
the photon exchange yields. 
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where Xi^'^^ is the tree amphtude for neutrino scattering off the assembly 
of free quarks, and Q = 1/6 is the mean charge of the up and down quarks. 
Note that the second term (proportional to Q) originates from the product 
of two antisymmetric tensors e'^'^'^f of leptonic and quark gamma matrices. 
Similarly for the Z-boson exchange, we obtain 

^<^"''-l|^(3<3 + ^)>o.Q)>«°'. (41) 

The first term is again due to the product of two antisymmetric tensors. 
The second term is all the rest. 

In order to connect the quark-level amplitudes with hadronic ones, we 
assume that the ratio of the tree and loop amplitudes for the neutrino quark 
scattering is the same as that for the neutrino nucleon scattering [7]. This 
is justified at least for the logarithmic divergent part of the correction for 
the Fermi transition by current conservation. In the next section we justify 
this for the logarithmic divergent part of the correction for the Gamow- 
Teller transition in the presence of partial conservation of the axial current. 
With this prescription we write the radiative correction to neutrino nucleon 
scattering at the short distance as, 



+ 3Qlog ( ^ 





We observe that the M dependence (upper cutoff) that appears in the long- 
distance radiative correction (37) is cancelled by the first term in the brack- 
ets in (42), which renders the ultraviolet divergence in the QED radiative 
correction finite [7]. Rigorously speaking, this universality applies only to 
the logarithmically divergent part, and the constant terms might receive 
extra contributions. Here we simply take the calculation with the point 
nucleon for the constant term. 

Note that the weak coupling constant is determined by the rate of muon 
decay, where the same box-type Feynman integrals appears. The formula 
used to determine Gp from muon lifetime 



Glml 8ml\ ^ 3m^ ^ 




25 o' 



(43) 



means that the effects of the box-type diagrams in muon decay is absorbed 
in Gp- Thus, when we calculate radiative corrections in terms of Gp, we 
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should subtract those included vn Gp'- the relevant contributions for muon 
decay are obtained by putting Q = — 1/2 in (42): 



m 



^ 3 
87r2 1 2 I M2 



w 




|A^(0)|2 

(44) 



Adding (42) to the QED corrections 5? and subtracting (44) , the Fermi 
part inner correction of (37) is 



0\r 




where the terms proportional to qa/ fv are collected in C^, 



3, /m2^ 



+ 



(46) 



As noted above, the M-dependence in the first logarithm in (45) is can- 
celled by the term ^log ^M2/m2^ in 5?^'. The other M dependent term 

(3g'A/2/y)log(M2/rn2) in SfJ that arises from the vcctor-axial-vcctor inter- 
ference, however, is hadron structure dependent, which is manifest by the 
fact that it is not cancelled by the electroweak counterpart 3Qlog(m|/M2): 
the coefficients of the two logarithmic terms agree only for an unrealistic 
value Q = gA/2fv- For the evaluation of these contributions, one separates 
long- and short-distance contributions [7, 9], leaving the M dependence that 
appear in (45) as the lower cut-off of the short-distance integral, whereas 
estimating the long distance contribution by rendering its logarithmic 
dependence milder with nucleon structure taken into account [9] . 

We can proceed in a parallel manner for the Gamow- Teller transition. 
The basic observation is the universality of the logarithmic divergence that 
also applies to the Gamow- Teller transition for the axial-axial contribution, 
as shown in the next section. With this universality we can identify the 
QED correction given in terms of nucleon matrix elements in (38) with 
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that of free quark theory, (40), and the match of the coefficients means 
the cancellation of the M dependences between the two expressions. We 
cannot prove the universality for the constant terms and those proportional 
to fv/SA- The corrections for the Gamow- Teller transition may have a 
more room of extra contributions from strong interaction than for the Fermi 
transition, unless chiral symmetry is exact. Again, we simply assume the 
calculation of the point nucleons as was done for the Fermi transition except 
for those proportional to fv IdA- We obtain for the Gamow- Teller transition. 



87r2 



(47) 



where the terms proportional to fv/dA collected in C*^^ are 



for point nucleons. 

For the vector-axial-vector interference term, we encounter exactly the 
same problem as for the Fermi transition. The two logarithmic terms, 
3Qlog(m|/M2) in (47) and (3/y/25A)log(MV"ip) in (48), match only when 
Q = fy jlgA, another unrealistic value. So, we take the same prescription 
carried out for the Fermi transition by Marciano and Sirlin [9], by includ- 
ing the long-distance part correction introducing nucleon form factors. We 
present the calculation of C^"^ in Sect. 6. 



5. Universality of the logarithmic divergent part for the 
Gamow- Teller transition 

Abers et al. [3] gave a proof that the logarithmic divergences '^ogiVP' jm^ 

that appear in the Fermi part, i.e., (5?^', is universal in that their coefficients 
are modcl-indcpcndcnt regarding the structure of hadrons: the coefficient 
of the logarithmic divergences does not depend on whether one deals with 
hadronic or quark current, in so far as these currents satisfy the generic 
commutation relations. Here, we argue that the same statement applies 
to the Gamow-Teller counterpart in b'^^ . The currents and the current 
commutation relations are given in Appendix D. We note that the proof 
described here is implicitly included in the work of Sirlin [16], which shows 
that the logarithmic divergences that appear in semi-leptonic processes are 
universal in the electroweak theory. We show explicitly the universality of 
the divergences between the four Fermi theory and the electroweak theory 
for the axial-vector correction to the Gamow-Teller part. 
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We introduce two types of Green's functions: 

Tx,.{k,VuV2) = i J d^xe^'-^pilT (t,(0)j---(x)) \p2), (49) 

Tx;,Ak,q,Pi,P2) = J d'xe^i-^J d'^ye^'^-y {p,\T {tx{x)j;"'- {y)jt-^- (0)) [p^) 

-STx^^u , (50) 

where J^'™ (a;) is the electromagnetic current and tx{x) is the V — A weak 
current 

tx{x) = Vx{x) - Ax{x). (51) 

The term —5Tx^v in (50) subtracts the pole term at P2 = rrip in Tx^v{k,p\,p2)-, 
which corresponds to mass renormalization. The considerations in what fol- 
lows does not depend whether the weak current (51) is expressed in terms 
of quark or nucleon fields, because our discussion is based solely on the 
current-current commutation relation between vector and axial-vecrtor cur- 
rents. 

With conservation of the electromagnetic current and current algebra 
applied to (49) and (50), we obtain identities: 

PTA^(fe,pi,p2) = {pi\txmP2), (52) 

k^Txf,ik,pi,p2) = {pi\t^iO)\p2)+Mi,-{p,-p2)^Tx^{k,puP2), 

(53) 

k^Tx^,,,{k.q,pi,p2) = -Txu{P2-pi-k-q,pi,p2)-k''STx^^, (54) 
q^Tx^v{k,q,pi,p2) = T^v{p2-pi-k-q,pi,p2) + T^,j.{k,pi,p2) 

+M^, - q^STx^,, (55) 

where 

= - J d'xe'^P'-P^-''y^p,\T{d ■ Aix)j;^-{0)\p2), (56) 

and 

M^. = -zj d^xe'^-^l d''ye'^-y{pi\T[d-A{x)f^'''-{y)jr-m) \p2) 

(57) 

arise in (53) and (55), as a result of non-conservation of the axial current. 
If we can set P2 — Pi ^ 0, then (53) is simplified as 

k^Tx^{k,pi,p2) = {pi\t^{0)\p2) + M^. (58) 
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By differentiating (55) with respect to and setting q = P2 — Pi, 
get the identity 



d 

T\ij,u{k,p2 - Pi,Pi,P2) + q^-Q^Tp,^u{k,q,pi,p2) 



q=P2—pi 



(59) 



q=P2—pi 



Neglecting again q = P2 — Pi, we write 



d 



d 

T\p.u{k,P2-pi,PuP2) = ■^Ti^u{-k,Pl,P2) + g^i^t^'^ ^ ^I'^^^Pt^^) 



q=P2-pi 



(60) 



V2 J (27r)^ 
Gv 2 f d^k 1 



Wc redefine A^(^) of Figure 1, including the effect of strong interaction 
in terms of the Green's functions (49) and (50): 

= -^e'/ ^^j^^.7'(l-7')5F(fc-^)7'^^^e(^)TA^(fc,Pi,P2) 

{J^\n{k,pi,p2) 

^Vy{k^-i^ + k^'-i^ - g^^'-i ■ k - is^f'^'-fal^kp - 2^7^^ (1 - ^^)ve{i) 

(61) 

after some manipulation of gamma matrices in the leptonic part. The first 
and second terms in the brackets of (61) can be expressed by using (52) and 
(58) in terms of the single-current matrix element. 

The third term requires a manipulation of g^^Txn{k,p\,p2)- It is known 
[12, 13] that the large energy limit of the Green's function (49) may be 
expressed in terms of the equal-time (ET) current commutator 

limfcO^+^r,^(A:,pi,p2) = ^ / d?xe-'k-^{p^\ [iA(0),j- --(^)l,^ IP2). (62) 



The current algebra relation, used in (62), 



ET 



iA(0),ir-(^) 



ET 



-2tQ{x)5^{x) 



(63) 



determines the high energy behaviour of g'^^Txn{k,pi,p2), which reads 



limfeO^+ooTA^(fe, pi , p2)g 



-^(pi|io(0)|p2). 



(64) 
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This is non-covariant looking. If we are interested only in the leading term 
in A;2, we can follow the trick of [3] to rewrite it in a covariant form: 

Vnnk^^+ooTxf.{k,pi,p2)g^'' - ^^ _ (pi|t^(0)|P2)- (65) 

The evaluation of the fourth term in (61) requires knowledge beyond 
generic current commutation relations. Writing explicitly the current in 
terms of nucleons or quarks (see Appendix D) , we obtain 



tx{^),jl'''-{x)\^^ = ^Q{kot%Q)-g"^k-m}5\x), (66) 



where Q is the mean electric charge of the fundamental doublet, i.e., Q = 
1/2 for proton and neutron, or +1/6 for quarks. The presence of Q means 
that this current commutator is not model-independent. Note also that the 
commutators used to derive this expression are those beyond ordinary cur- 
rent algebra; they are determined only with the aid of models of hadrons 
(see Appendix D). We see that the antisymmetric tensor in the l.h.s. of 
(66) converts the vector (axial-vector) currents into the axial-vector (vec- 
tor) currents. Therefore, the qa dependence appears for the vector current 
contributions upon taking the hadronic matrix element. This contrasts to 
the situation for the first three terms in (61), which are determined by 
current algebra and conservation laws; the effect of strong interaction for 
the axial-vector current is properly absorbed in after taking the matrix 
element with nucleons. 

We are interested only in the leading high energy behavior and again 
write (66) in a covariant form: 

-i limfc2^+<^e^'""^fepTA^(fc,pi,p2) 

- 40 (,--^) (,,1^.(0)1,,)^,-^. (67) 

The fifth term in the brackets of (61) does not contribute to its divergent 
part. Thus, the divergent part of (61) is evaluated as. 



divergent 



+4Q U-^)iP^\tM\P2)^^. 



P2 
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Gv 2 

~ 

V2 



1 



2 + 



1 



(27r)4 fc2 - A2 (fc - If 
(piKa(0)|P2) + Ma 
A:^ 1 



mi 



+3Q(pi|iA(0)|p2) 



167r2 



logM^ + 3gA^ 



(0) 



where 



A^(0) = ^V.^\l - 7')Veie){piM0)\P2) 



(68) 
(69) 

(70) 



is a generalisation of (3) including the strong interaction effect. We have im- 
plicitly assumed above that Mx in (68) does not contribute to the divergent 
part, because axial current conservation is broken only by soft operators. 

The external nucleon-line renormalization also receives strong interac- 
tion effects. Using the same notation A^^*) as in (26), while including strong 
interaction effects, A^(*) is expressed in terms of the Green's function (50), 



Gve^ 



2^2 



(27r)4 A;2 - A2 



(71) 



The last term is external line renormalization of the lepton. 

The integral on the r.h.s. is expressed in terms of T^i, by using (60) and 
further (pi|i^(0)|p2)- Retaining only the divergent part, and using (65), we 
find 



2^2 



111/ 



(27r)4 fc2 _ X2 

d d 



dk^ 



q=P2-pi 



' Gve' 



d^k 



2^2 

d 



J W 



(27r)4 P - A2 
2ki' 



dk^ \B + 2k-p2 



V,^^{1 - 7')i;e(^) 



d 



(Pi|i/.(0)|P2) +5^'^^(AV - q'^Tp^u] 



q=P2-pi 



(72) 
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Again, axial-current conservation is broken only by soft operators and the 
M^i, term does not contribute to the divergent part of M^^"^ . Thus we can 
drop M^y — qfSTpniy. We, therefore, retain only the first term of (72) to 
evaluate the divergent part of Al^*); 



divergent IGtT 



„2 



logM^TWW + jy^K) - l| 



rlogM^TW^^l (73) 



167r2 

Adding (69) to (73), the divergent part of M^'"^ + M^^^ becomes 

/3 



,. ^ 167r2 

divcgcnt 



logM^ + 3g7W(o)^ . (74) 



The first part in the brackets of (74) gives divergent terms that are pro- 
portional to fy and g\ in the amplitude squared, as determined by current 
algebra, conservation of the vector current, and softly-broken axial current 
conservation. Hence, we can conclude that the coefficient of the logarithmic 
divergence ^log{M'^ /nip) in (37) and that in (38) are both independent of 
the model of strong interactions, and smoothly continue to the electroweak 
calculation. On the ohcr hand, the second term of logarithmic divergence 
which arises from the interference of the vector and axial-vector currents 
are model-dependent; the appearance of Q is a manifestation of the model 
dependence. 



6. Long-distance contributions to the vector-axial-vector 

interference terms 

To evaluate the long-distance contribution to the axial-current correction 
to the Fermi transition, Marciano and Sirlin [9] proposed the prescription to 
make the logarithmic 'divergence' milder by introdticing nucleon form fac- 
tors in the evaluation of vertex corrections. The result does not contain the 
cut-off mass M any longer. The calculation for the Gamow- Teller transition 
is carried out parallel with that for the Fermi transition. 

Marciano and Sirlin [9] considered only the vector and axial-vector cur- 
rent contributions, implicitly assuming that the contribution from weak 
magnetism is suppressed by the proton mass 0(l/mp). We find that, since 
the loop momentum becomes of the order of the form factor mass, typically 
nip (p meson mass), an 0{nip/nip) contribution may be non- negligible. This 
contribution turns out to be particularly large for the Gamow-Teller tran- 
sition. 



22 



TEXT PRINTED ON FEBRUARY 1, 2008 



We work with the effective electromagnetic vertex, 

j''Fi^\k'') - J-a^'-k,F^\k''), (75) 

for the proton, where Fi^\k'^) and F2^\k'^) are the form factors, and 

similarly for the neutron with the form factors, F^^\k'^) and F2"'\k'^). 
The form factors are also included at the weak vertex of the nucleons, 
Fy {{p2 — k — pi)'^)) ~ Fv{k'^) for the vector vertex and Fa{{p2 — k — pi)"^)) ~ 
FA{k'^) for the axial- vector vertex. We also retain in the weak vertex the 
weak magnetism term with the form factor Fy/ik"^) via 

" -(rxp{P2 -k- pi)PFw{{p2 -k- p^f « -^axpk^Fwik^), (76) 



2mN 2mjv 
where m]\f = {mp + mn)/2. The effective nucleon vertex Ri^\{pi,p2,k) 
defined in (20) that contributes to M^'"^^ then reads, 

Rf^xiPi,P2,k) = ^{axpkP){1p2-k)^F[^^{k'')Fw{k^) 

+ {-m [fvFvik') - gAl'FAik^)) + 2^i^>^p'''')Pw{k^)] 
xia,.kn{Fi'\k') + F^\e)}. (77) 
We now decompose the inner part as, 

where 

d^k 1 11 



" ^/2 J (27r 



V2 J (27r)4 (A; - - m2 (p2 - A;)2 - m2 fc2 - A2 
Xi;^7^(l - 7^) {{k - 2iY + ia'^'k,,} Ve{i) 
xuniPihx {fvFvik^) - gAl^'FAik^)} {a^pkf)up{p2) 

x{Fi^\e)+F^''\k')}, (79) 
d^k 1 11 



(27r)4 (A; - - ml {p2 - kf -mjk'^- 
xv,j\l - 7^) {{k - 2£r + ia^^'k,} Ve{£) 

{axpkP){2p2-k)^F[^\k') 



-ii'^xpknia.ukn {Fi''\k^) + Fip\e)} 



Up{p2)Fwik^) (80) 
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and likewise for the neutron, M.n^'^^^ and A^i"^''"™''. Eq. (79) is the 
term due to the V and A current interaction and (80) is due to the weak 
magnetism. The calculation of the four terms in A^("^) of (78) is given in 
Appendix E. 

Results for TWj,"^'^^^ and M^n^'^^^ 

We are interested only in the /y^^i-terms and we subtract the pure 
vector and pure axial- vector parts that are proportional fy and g\. We 
find 



spin 



fv9A 



= 32GvmnmpEEi, I ^ j fvgA 



6 (C^P'^^ + 4"'^)) (1 + /?cos^)(l)^ 



where and C-r are defined by the integrals 



(81) 



167r2 



mt 



(82) 



kn ky 



(27r)4 (P)2{(p2 - fc)2 



167r2 



-^{F,(^)(fc2)+Fi^)(A:2)}F^(A;^) 

<7;..C(f + (P2)M(P2).C(f I , (83) 



and arc evaluated numerically using the empirical dipole nucleon form fac- 
tors (see Appendix F) as, 



Cipy) = (1 + /xp) X 0.240 = 0.671, 
Cipy) = (1 + /xp) X (-0.198) = -0.554, 
Q{p,A) = (1 + ^p) X 0.261 = 0.729, 
C(P,A) = (1 + pip) X (-0.212) = -0.592, 



(84) 
(85) 
(86) 
(87) 



where fip = 1.793 and /x„ = —1.913 are the anomalous magnetic moments 
of the proton and the neutron. 
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The evaluation of the neutron's counterpart C'J'^ , C^*'^, CJ'^ and C^'^ 
goes similarly, with the result, 

C^^y) = /i„ X 0.240 = -0.459, (88) 

4".^) = /i„ X (-0.198) = 0.379, (89) 

^ ^ Q 261 = -0.499, (90) 

4"'^) = n^x (-0.212) = 0.405. (91) 

The first terms in the square brackets of (81) are the axial-vector current 
contribution to the Fermi transition, and the second terms are the vector 
current contribution to the Gamow- Teller transitions. The computation of 
the former terms (i.e., those proportional to (1)^) was already done by the 
authors of [9] and [10], and our evaluation 

fv9A X 6 (ciP'^'^ + C("'^)) = 2 X 0.875 = 1.75 (92) 

agrees with their results allowing for the difference that arises from the 
different values of the input data ^. 

Let us remark that (84) - (87) and (88) - (91) reduce to 



C^,v),C^P,A)^U,J!^\+l (93) 




C^^),^"'-^) ^0, (94) 
C(py)^C^A)^_^^ (95) 

C{n,V)^C{n,A) (96) 

for the point-like nucleon with a vanishing anomalous magnetic moment, 
and (81) reduces to the /y^^-terms in (22). 

Results forM^f^""""^ and M^f^""""^ 

The calculation is also straightforward for the effects of weak magnetism, 
though it is rather lengthy. We refer the readers who are interested in the 
calculation to Appendix E; here we only record the final result: 

spin 



^ Our value 0.875 in (92) should be compared with 0.885 of Marciano and Sirlin [9] 
and with 0.881 of Towner [10]. The slight difference is due to different choices of qa 
and ruA in the form factor. 
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ruN 



32G'^innmpEE^ 
+gA {rUpS^P^ - m„f(")) (3 - pcose)^{a)' 
+lfv {mp£^^ + m„f ('^)) (1 + /3cose)(l)2 



where the constants and are defined by 

d^k k^'k•' 1 



2gA (mpP^f) + m„p(.")) (3 - pcosO) ha)' 



(97) 



167r2 



(98) 



167r2 



P2 



(99) 



and are evaluated, by using the phenomenological form factors of the proton 
in Appendix F, as 



= _ X 0.240 = 0.890, 
= (|Xp - /i„) X (-0.198) = -0.735, 
£(P) = (/ip-/i„)(l + /xp) X 0.0836 = 0.866. 



(100) 
(101) 
(102) 



[Note that vi^^ does not appear in (97).] We also define the corresponding 
quantities for the neutron and obtain an expression similar to (97), where 



= 0, 

f H = (^p - //„)//„ X 0.0835 = -0.592. 



(103) 
(104) 
(105) 



Summary of long-distance contributions 

The long-distance contributions of M.^'"^^ defined by (78) is summarized 

as 

spin 
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32GlmnmpEE^ ■ ^ 



C^il + Pcos9)f^{l)' 



(106) 



Here the constants and C^""", which were introduced in (45) and (47), 
are obtained from (81) and (97): 

f^C^ = GfvgA (C^^) +C(--^)) + ^ (m,^(^) +m„f H) , (107) 



2mj\[ 

, 9 A 



ny) 



rriN 



{2 (mpP^) + m„Pi")) + {mp£^^ - m„f } , (108) 



which are evaluated numerically as, 



= 1.751 + 0.409 = 2.160, 



C^"^ = 0.727 + 2.554 = 3.281, 



(109) 



where the two parts represent contributions from the (V,A) interaction 
and weak magnetism. Note that the weak magnetism contribution is non- 
negligible for C^, and is even larger than the (y,A) contributions in C'^'^. 



7. Summary 

Full one-loop radiative corrections are calculated for neutrino-nucleon 
quasi-elastic scattering for both Fermi and Gamow- Teller transitions. We 
separate the corrections into the outer and inner parts a la Sirlin. The outer 
part is infrared and ultraviolet finite, and involves the positron velocity. This 
part takes different forms for angular independent and dependent parts of 
the differential cross section. The calculation of the inner part requires a 
scrutiny regarding the continuation of the long-distance hadronic calculation 
to the short-distance quark treatment and the dependence on the model of 
hadron structure. We have shown that the logarithmically divergent parts 
do not depend on the structure of hadrons not only for the Fermi part, 
but also for the Gamow- Teller part. This observation has enabled us to 
deal with the inner part for the Gamow- Teller transition in a way parallel 
to that for the Fermi transition. The inner corrections contribute to the 
angular-independent and dependent parts of the cross section in the same 
manner, so that they are written as universal multiplicative factors on the 
coupling constants, fv and gA [4]. 
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The resulting radiative correction to the differential cross section (5) 
takes the factorised form: 

A{P) = {l + 5,^,{E)}(f^{lf + -g\{a)^) , (110) 
B{P) = {l + 5out(^)} - lfA{<rf) , (111) 

where 5out(-£') and 5out{E) are given by (36) and (39), respectively, and 
the inner corrections are absorbed into the modification of the coupling 
constants, 

fj = /^(l + C), 

gl = 9a{^ + S?J) , (112) 

with Sf^ and 6^ given by 



^GT _ « \^^^J'^z\ , .(mp\ , , , ^gt 



and the constants C^'*^""^ by (109) [a is the fine structure constant, a = 
e^/ (47r)]. If we set the mass scale, M, at which the short- distance asymptotic 
behaviour onsets, to be 1 GeV, we obtain 

5l = 0.02370 - 1.16 X IQ-^ x log (j^) 

= 0.02370 ± 0.0008, 

5^^ = 0.02616 - 1.16 X 10"^ x log (^^y) 

= 0.02616 ± 0.0008, (115) 

the error corresponding to a multiplicative factor of 2 in the scale M. 

We present numerical values of 6out {E) and 5out {E) in Table 1 and Figure 
3, where T = E — rUe is the kinetic energy of the positron. In particular, at 
the threshold E = rUe, they are given by 

f 27 3 / ^n'^ \ 1 
<^out(me) = + 2^og = 0.01835, (116) 



^out(me) = + llog I = 0.02067, (117) 
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and for a large positron energy, they behave as 

5out{E) and 6ont{E) differ only by « 10%, and and also differ by 
10%; hence the radiative correction to the angular dependence is at a 0.1% 
level. 



T (MeV) 


/5 


'^out(^) 


5out(^) 


0.01 


0.1950 


0.0180 


0.0203 


0.05 


0.4127 


0.0173 


0.0193 


0.1 


0.5482 


0.0166 


0.0184 


0.2 


0.6953 


0.0156 


0.0171 


0.3 


0.7765 


0.0149 


0.0160 


0.4 


0.8279 


0.0142 


0.0152 


0.5 


0.8629 


0.0137 


0.0145 


0.6 


0.8879 


0.0132 


0.0139 


0.7 


0.9066 


0.0129 


0.0134 


0.8 


0.9209 


0.0125 


0.0129 


0.9 


0.9321 


0.0122 


0.0125 


1.0 


0.9411 


0.0119 


0.0121 


2.0 


0.9791 


0.0099 


0.0096 


3.0 


0.9894 


0.0086 


0.0081 


4.0 


0.9936 


0.0077 


0.0070 


5.0 


0.9957 


0.0070 


0.0062 


6.0 


0.9969 


0.0064 


0.0055 


7.0 


0.9977 


0.0058 


0.0500 


8.0 


0.9982 


0.0054 


0.0045 


9.0 


0.9986 


0.0050 


0.0041 


10.0 


0.9988 


0.0047 


0.0037 


20.0 


0.9997 


0.0023 


0.0013 


30.0 


0.9999 


0.0009 


-0.0002 


40.0 


0.9999 


-0.0001 


-0.0012 



Table 1: Numerical values of the outer corrections, Sout{E) and 6out{E) as 
a function of the kinetic energy T = E — rrie of the recoil positron. 
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Fig. 3. The outer corrections, ^out(-£') (solid line) and 5ont{E) (dotted line) as a 
function of the kinetic energy T = E — rrie of the recoil positron. 



The dominant part of the inner correction arises from 41og(m^/mp), 
which amounts to 18.18. The change of the scale M = I GeV in log(mp/Af ) 
by a factor of 2 leads to an error of ±0.69. The model-dependent, long- 
distance contribution from axial- vector- vector interference terms contributes 
by 2.160 for the Fermi transition, and 3.281 for the Gamow- Teller transition 
(there is an additional contribution of unity for the Gamow-Teller transi- 
tion). We find a sizable contribution to the axial- vector- vector interference 
terms from the weak magnetism, which has been ignored in [9] and [10]. 
This contribution is even larger than that of the {V, A) interaction for the 
correction to the Gamow-Teller transition. If we would use the point nu- 
cleon calculation for the Gamow-Teller transition as an extreme case the 
long-distance contribution to C*^^ will be 1.1. We may, therefore, conclude 
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that the uncertainty from the model dependence is conservatively no more 
than 10% for Si^. 

The formulae we derived are applicable not only to neutrino-nucleon 
quasi elastic scattering, but also to other neutrino reactions off nuclei. 

The radiative correction to nuclear beta decay is also described by (110) 
with (5out replaced by the well-known g function of [1, 4] 



^out 



a 
2^ 



3 In 



1 



nip 



4 + /3^U + /3 



+4 ( -tanh-^^ - 1 
+^tanh-i/3|2(l + /32) + 



Eq — E 



3E 

(Eq - Ef 

6£;2 



^+ln^ 

2 vfip 



E) 



4tanh-^/3 



, (120) 



where Eq is the end point energy of the electron. The difference from the 
case of 0+ — *■ 0+ decay is that we have now somewhat larger inner part 
correction for the coupling constant that is multiplied on the Gamow- Teller 
matrix element. For example, the neutron decay width is given by 



r = 



(121) 



where /{Eq) = 1.71483 [18] includes the outer radiative correction. The 
value in the square bracket is 5.966, which is 0.20% larger than 5.954 that 
would be obtained by assuming 6^ would equal [19], where gA = 1.2670 
[20] is used as a provisional value. Unfortunately, the values of gA reported 
in the literature or in the table of Particle Data Group [20] are neither the 
bare value nor the value including the outer radiative corrections defined 
through the prescription discussed in this paper. Therefore, the current 
accuracy of gA does not warrant the use of our formulae, e.g., to determine 
precisely \Uud\ = cosOc'^, and we should postpone a detailed numerical 
analysis until an accurate estimate of gA will become available. 



^ Particle Data Group [20] gives a value gA = 1.2670 ± 0.0030. This value is obtained 
by naively averaging over values reported in the literature, one of which is obtained 
including the outer radiative correction (up to an unkown constant; the inner correc- 
tion is ignored) , and some others are results which do not include radiative corrections 
at all. Unfortunately, the value with radiative correction reported by [22] differs from 
others that ignore radiative corrections [21] by an amount much larger than the dif- 
ference that would be accounted for with radiative corrections, and hence much of 
this difference is probably ascribed to other systematic errors of experiments. For 
this reason we give up to apply corrections to find a bare qa from the available data. 
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Our final remark is that one may not necessarily need the formulae we 
derived for the purpose to obtain a precise estimate of neutrino nucleon 
scattering cross section or of |f7ud|) if 9a would be properly obtained by 
using the outer-part radiative correction formula in which the inner and 
outer corrections are accurately separated and the inner corrections are 
included in qa- The true power of our formulae emerges when one would 
deal radiative corrections for processes involving the neutral current. 
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Appendix A 

Calculation of M^'''^^ in (17) 

We sketch the calculation of M.^'"'^^ in (17). The integration formula 
necessary for (17) is given in Appendix of Ref. [3], 



(27r)4 {k - £y -ml(k- p^f - mlk'^ - 

1 1 / rr? ^ 

--tanh-i/3 + -logl ^ 



1 1 



+ + ^)tanh-i/3 - llog j I . (A.l) 



Putting (A.l) into (17) and taking the nucleon static limit, we get 

J^tanh-/? 

x[v,j\l-j^)ia'"'ve{i)KK{Pi)Wx{pi,P2)up{p2)] , (A.2) 



where only the term proportional to nip in (A.l) is retained. The spin 
summation of the tree-one-loop interference amplitude is expressed as 



(0) 



spin 
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(A.3) 



where 



spin 

Q^P = J2[v,j\l - j')ia''''vemAv.Y{l - l')ve{i)r, (A.4) 

spin 

are evaluated as 

Kqo = SfyUinmp, Kij = 8g\mninpSij, (A.5) 
qOO = 8E'^E^P{(3 + cos9), Q'^Sij = 8E^E^P{3p - cos9). (A.6) 

We arrive at 



spin 



= 32G'^mnmpEE^ [t^] tanh'V {fviP + cosO) + g\{3P - cos6')} , 



(A.7) 



which is the M^^^) p^rt of (18). 

Appendix B 

Calculation of M^'"^^ in (19) 
Neglecting A and rrig, Ai^""^^ is written as 
Gv 2 f d'^k 1 1 



V2 J (27r)4 (fc2)2 - ky - ml 

XUnipihxifv - gAl^)<ytipk^Up{p2) 



Gyf^ 

V2 Vl67r2 



T PX" 



,7^(1 - ^^)ia''^Kve{t) 



(B.l) 



where 
i 



167r2 



' " - J (27r 



(Pk 1 



(27r)4 (A;2)2 (^^ - kf - ml 



167r2 



1 /m2\ 3' 



(B.2) 
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and 

X% = - l^)<^'"'Ve{mUn{pihx{fv - 9Al^)a^,pUj,{p2)]. (B.3) 

We shuffle ^\(Tfxp in (B.3) using 

^^YY = 9^^Y + g^""!^ - g^^J^ + ie^'^^^Jpl^ (B.4) 

in such a way that a single gamma matrix is sandwiched between spinors; 
we find in the nucleon static limit 

/"X.p = -6{fv + gAnVul^'il - l')Ve{mMPlhoMP2)] 

(B.5) 

l^{p^y{p^)PX,p = -3fv\v,-''' 
rUp 

+ (2/y + 9a)[v.Y{1 - l^)Ve{imn{piha^Up{p2)]. 

(B.6) 

The first and second terms correspond to the Fermi and Gamow- Teller parts, 
respectively, in (B.5) and (B.6). The angular distribution is determined by 
the leptonic parts, 

spin 

Y,[v,j\l - j^)vemyul'{l - l^HiiMj = 8EE,{3 - /3cos/3).(B.8) 

spin 

Combining with the nucleon parts, 

^[Un{Plh'^Up{p2)][Un{pihoUp{p2)]* = SmnUlp, (B.9) 
spin 

^[Un{Plhi7^Up{p2)][UniPlhj7^Up{p2)]* = Smnnipgij , (B.IO) 
spin 

we end up with the A4 ^""^^ contribution to the differential cross section given 
by, 



spin 

= 32G^ymnmpEE^ f ^ J (i + pcosp) 
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6/y ifv + 9a) { ^log 



M2 



mt, 



+ 77 +3/ 



+Z2GvmnmpEE^ 



87r2 



(3 - /3COS/3) 



65A [fv + 5a) l^log + ^1 + 5a(2/v + gA) ■ {-\ 

(B.ll) 

This is (22). The first and second terms correspond to the Fermi and 
Gamow- Teller contributions, respectively [the matrix elements are explicitly 
denoted in (22)]. 



Appendix C 

Calculation of bremsstrahlung 

The bremsstrahlung diagrams contain infrared divergences and the Feyn- 
man integrations require careful treatments, in particular, to calculate the 
angular dependent part. The terms in the first line of (29) do not contain 

and are angular-independent; hence, the angular integration of k is easy. 
Those in the second line of (29) are linear in and are proportional to 
cos^. After some manipulations, (28) becomes 



da{ve + p — > e"*" + n + 7) 



d{cos0) 
47r2 



TT 



(/y + 35i)|-^'tanh-i/3-log( 



A 



E — rUf 



6 

k=l 



k=l 



(C.l) 



Here the integrals Zk, which do not contain cos^, are given by 



J nip 



E~E' ^\E{1 + P)J 

2{E + me 



E' + \£'\\ 



(tanh-^^)^! , 



(C.2) 



pE-\ 

I2 = -E dE'log 



rup 
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= E"^ (-tanh- V + 13) , (C.3) 

2 Jme \ rUe J 

= - ^] , (C.4) 



li = - 



e; fE-X 1 



5 



Xr = -m^£; / dE' 



2 ^ Jm, \2E'{E-E>) + 2\k\\i\+\^ 

1 



2E'{E - E')-2\k\\l\ + X^ 



fE-\ 1 /2£;'(S-^') + 2|fc||^| + A2 



me 



""l 2^'(^-^')-2|fc||£| + A2 E' + 



{^L - ^(tanh-i/3)2 + tanh-^/? ■ log2| . (C.7) 

Collecting (C.2) to (C.7) together, we obtain the angular independent part 
(31). 

The integrals in (29) that are proportional to cos 9 are given by 



= -!-/ dE'\e\r-^5{E-E'-u) 



A2 



{2k-l' + X^Y \ a;2 



= E^^cos^jl + ^ - ^^^^^ - (1 + /3)tanh-V 



(2A; • i' + A2)2 r a;2 



J,^^ dE'\e'\ / — 5(E - - a; 

27r£;,, Jme ' ' J LO ^ 
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E'^Pcosel^tanh.-^P -(3+ ^^^(tanh-i/?)4, 



(C.9) 



J3 



= _1_ / dE'\£'\ / —S{E -E'-io) 
= E^/3cose|^tanh-V 



{2k ■ £' + A2)2 



,-i.,2 13/3 vlHZ 



(C.IO) 



J4 = 



2nE^ 



E-\ 



dE'\l' 



f 1 ^ - fP' 

= E^pcose I -tanh-1/3 ^{tanh-^pf - 



1 i-E-X r d'^h 

J, ^ dE'\t'\ r^5{E-E' 



CO 



I , (c.ii) 

{2k ■ I' + A2)2 a;2 



= E^l3cose I -^tanh-i/3 - ^(tanh" V)^ + ¥ } • 
I 4 8/j 8 I 



(C.12) 



Collecting these terms, we get (32). 

Let us now explain the treatment of infrared divergences that appear in 
the integral of J7fc, taking J\ as an example. After angular integration of 
the k variable, the integral is 



Jl = COS0 



E-X 



dE'\e'\ 
1 



2 \2E'{E - E') - 2\k\\i'\ + 

} 



2E'{E-E') + 2\k\\e'\ + X'^ 

\k\\e'\ 

'2{E-E'f 

E' (2E'{E-E')-2\k\\e\ + \^^ 

'2{E-E') °^ \2E'{E-E') + 2\k\\£'\ + X'^^ 



(C.13) 



A special care is needed because there are two sources of the photon mass 
(A) dependence. One comes from the the integrand and the other from the 
integration region. To handle this problem, we split the integral into two 
terms, as Ji = J[ + {J\ — Ji), where 



Ji = cose / dE'\£' 

J me. 



E' ( E'-\e'\ 

E-E' 2{E - E>) °^ I + l^'l 



(C.14) 
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which is defined by putting A = in the integrand of (C.13), while retaining 
A that appears in the integration region of E' . (Recall the relation = 
^E''^ - mi and |fc| = ^{E-E'f - A^) . 

Integral (C.14) is readily handled by changing the integration variable 

as 



1 



(C.15) 



with the result, 

J[ = E'^l3cose 



3 2 

4^ /3 /3 



1 



/3 



-+15] tanh-i/3 +[- tL+ ] (tanh-i/?)^ 



4/3 



2/3 



+ 2L 1 



+ (/? - tanh" V) log 



'1-/3' 
1 + /3, 



me 



1 + 



(C.16) 



The evaluation of Ji — J{ requires a great care, since its integrand is 
non- vanishing only in the vicinity of the edge point E' = E — \. To perform 
this integration, we choose the variable, 



7/ 



E- E' + ^'\E - E'Y - A^ 



(C.17) 



In the limit A — > 0, 

Ji-Ji = E^(3cose 



^{tanh-^Pf + ^ + ^tanh-i/3 



(/3 - tanh-i/3) log2 + h (- 



2f3 



+ /9 



(C.18) 



which is infrared finite. By adding (C.18) and (C.16), the integration of jTi 
is completed to give (C.8). 



Appendix D 

Currents and their commutation relations 
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We present explicitly the currents and the current commutation relations 
that appear in Sect. 5. Consider a fermion field ijj, which belongs to some 
representation of the internal symmetry group and satisfies the equal time 
(ET) canonical commutation relations ^/^^(y) r ~ ^af^S^i^ — y)- 

L ' J ET 

The vector and axial-vector currents are constructed as, 

j^ix) = ^{xh^T-i;{x), jf{x) = ^(x)7M7'r»V(x), (D.l) 

where T'^ are the generators of the internal symmetry group. 
By applying the canonical commutation relations, we obtain 

-ij{y)lul^l^.T''T''^jj{x)5\x - y), (D.2) 

VS(a^)7^7'7°7..r»rV(j/)<5'(a; - y) 

-V^(y)7.7%M7'r^rXx)(5^(a; - y), (D.3) 

Vi(a;)7/.7Vr«rV(y)53(a; - y) 

-i^{yh,-f\T^T''^{x)S^{x - y). (D.4) 

Prom (D.2)-(D.4), conventional current algebra follows for v = Q. 

The commutation relations we use are derived from (D.2)-(D.4). One 
class of the commutation relations is 



ET 



ET 



ET 





ET 






ET 




Jl\x),jl\y) 


ET 



-2V'(a::)7o 



Tp{x)5^{x - y), 



= -2^{x)-foj' [t«, T"] ^{x)S^{x - y), 
= -24^{x)jo [t%T^] 4^{x)S^{x - y). 
For the other class of the commutation relations, by noting 



£^'^'^''7m707. 



2i {gh'a 



9, 



11 



we obtain 



'j;ix),jt{y) 

= 2i [gU - 9a9'o) ^{xh'^l' {r«,r^} ^{x)S^ix - y), 



ET 



(D.5) 
(D.6) 
(D.7) 

(D.8) 
(D.9) 



= 2^ [9^9'a 

'f;{x),jl\y) 



ET 



Vi(x)7'^{r^r^}V(x),53(a;-y), (D.IO) 



ET 



= 2i (g^gP, - g^gP) ^{x)^^' {r",r^} ^|;{x)^^{x - y). (D.ll) 
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Let us consider the case of SU{2) and suppose that V belongs to the 
doublet representation. The electromagnetic and weak currents are 

= ^jf.T'^iP, tx = V^TaT^V - V^7A7'r-V', (D.12) 

where the charge matrix and T" are defined by 

Q+ \ ^- _ / 



the electric charges of upper and lower components of the doublet are de- 
noted by and Q_. 

The commutation and anticommutation relations between and T~ 

are 

T", T^] = (Q+ - Q_)r- = T-, (D.14) 
{t-,T^} = {Q+ + Q_)T~ = 2QT-. (D.15) 
The relations arc translated as 

9^"" [tA(x), jTMet = -2*o(rc)<53(a; - y), (D.16) 
s^'^''''[t^{x),jr{y)]ET = -^iQ{9U- 9l9^)t%x)5\x-y). 

(D.17) 

The relation — (5_ = 1 holds independently of the model of the con- 
stituents of hadrons; hence (D.16) is model- independent. On the contrary, 
the mean charge Q depends on the constituents, e.g., whether ip denotes the 
nucleon- or quark-doublet (0 = 1/2 or 1/6), so (D.17) is model-dependent. 

Some commutators of the local currents may in principle receive model- 
dependent Schwinger terms. Within the standard model, however, the 
Schwinger term is c-numbers, and it does not contribute to beta decay of 
nucleons, as argued by Sirlin [8]. 

Appendix E 

Calculation o/>i(^3'^^) and >l(^3,«;m) 

To shuffle the gamma matrices in (79), we use 
\ {ta {fvFvik') - 9AFA{k^h') 

= -sx^pa^l^ {fvFvik^) - 9AFA{k^h^) , (E.l) 



7A (fvFvik^) - gAFAik^h") ,cT^p\ 

= i (5Am7p - 5Ap7m) [fvFv{k') - gAFAik'^h^) , (E.2) 
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and simplify nucleoli's gamma matrices in (79) using 

+1 [ta (fvFvik^) - gAFAik^h') , a^p] ■ (E.3) 

Applying these formulae and using the definition of Ca,T in (82) and (83), 
(79) is re-expressed as 

vrip 

(E.4) 



I 

nip 



We are interested only in those terms that are proportional to fv9A 
in the amplitude squared. By inspecting the gamma matrices in (E.4) in 
the nucleon static limit, we see that only the first and second terms in the 
brackets of (E.4) are potential sources of the fv9A terms; the third and 
fourth terms give those proportional to fy or 9\. The spin summation 
relevant to the first and second terms are 

E ex,ua[Un{pi)ri'' (/yC(f'^) - 9Al'Cj^''^^) Up{p2)] 
[Un{Plhp (fv - 9A7^) Up{p2)]* 
spin 

38AmnmpEE^fv9A {c^'^\l + PcosO) + C^'^\3 - pcosO)} , (E.5) 

E -^^X,raiP2UP2r[Un{Plh''l' (/y4^'''^ " ^At'^^'^^ ^^^2)] 

mp 

[Un{Plhp (fv - 9A1^) Up{p2)]* 



spin 

X 



spin P 
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spin 

12%mnmpEE^fvgAC^'^\'3> - PcosO). 



(E.6) 



It is obvious that (l + /3cos^) and (3 — /3cos6') terms correspond to the Fermi 
and Gamow- Teller parts, respectively. 

The interference with M.^^^ turns out to be 



J2 (a^(-3.^^)a^(o)* +A^(-3,VA)*^(0)^ 



spin 



Z2GvmnmpEE^ I ^ 1 fvQA 



-(6C(.f'^) + 24^'^)) (3 - l5cose)\{a) 



fv9A 



(E.7) 



where only the fv9A terms are retained. The same procedure applies to the 
neutron amplitude, and we arrive at (81). 

For the weak magnetism terms, the gamma matrices in the last line of 
(80) are simplified by applying 



{(Txpkni(T^,.n = \{axpkP,a^,k^} + ^[ax,kP,a^,k''\, (E.8) 

[axpk'', a^j.k"] = -ikxa^pk^ + ik^.axpk'' - ik^a^,^ + iex^uKk''{l ■ k)-f'''j^, 

(E.9) 

{c^XpkP, a^,k^} = -2 (kxk^, - k^g^^x) ■ (E.IO) 

Note that (E.9) and (E.IO) are antisymmetric and symmetric under the 
interchange of the indices A and fi. We also simplify the gamma matrices 
in the leptonic sector in (80) by using 



\l--f^)a^'''k^ = 1^ (A;V + 'tV) -5^^(7-^)|(l-7^) 

+ 1^ {k^^^ - - ie^'""'7pl^k,^ (1 - 7^). 

(E.ll) 

The first (second) line is symmetric (antisymmetric) under the interchange 
of A and fi. Neglecting rUe, I and A, and using (E.9), (E.IO) and (E.ll), we 
simplify (80) as 
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= -^Gve" ( I 7^77^7 ^ .F'f\e)Fw{k'') 

v/2 \2mNj J (2^)4 (fc2)2 - A:)2 - m2 ^ ^ ' ^ ' 

+2[V^J^{1 - 7^) {k-p2 + i(J^"'k^{p2)^) Ve{e)\[Un{pi)(TXpkPUp{p2)]^ 

^ 2( i \ f d'^k 1 1 

^\/2^'' V2m^ j V 2 y y (27r)4 (P)2 - fc)^ - 

x{F^\k') + F^\k')}Fw{k') 

x| - 6[t;^7 • fc(l - l^)Ve{(.)][Un{pi)Up{p2)\ 

+2e^>'^fk,[v^.^p{l - l'')vem^n{pi)(rx^up{p2)]^. (E.12) 

In (E.12) the term containing [un{p\)(T\pk'^Up{p2)]{k'^ — 2k ■ P2) does 
not survive the symmetric integration over k. In the static nucleon limit, 
we retain only those terms containing Un{p\)up{p2) for the Fermi part and 
UniPi)o'ijUp{p2) for the Gamow-Teller part. Using the definitions in (98) 
and (99), (E.12) is rewritten 



Gv / TUp 



\/2 167r2 yrriN 



+ %T#2 f — V^'M ?[^-7°(l-7'K(^)][nn(pi)np(p2)] 



■v/2 167r2 \m]\f 

-\ E e°'''[^.7fe(l-7')^e(^)][«n(pi)^i,-Up(p2)]|. (E.13) 

i,j,k=l ) 

The following spin summations are employed to calculate the interfer- 
ence with TVlW: 



^[UniPl)crijUpiP2)][Un{Pl)Y{fv " 9Al^)Upip2)y 



spin 

X 
spin 

= l2%gAmnmpEEy{2, - PcosO), (E.14) 
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^[Un{pi)Up{p2)][Un{pi)Y{fv - 9 A't')Up{p2)]* 



spin 



spin 



= GifviTinmpEEu^l + (3cos9), 



(E.15) 



spin 

X '^[Un(.Pl)crijUpip2)][Un{pi)Y{fv " 5A7^)«p(P2)]* 
spin 



= —128gAmnmpEEi,{3 — (icosO). 
We obtain for M.'^p^''^'^^ 



(E.16) 



spin 

( 1 
= 'i2GymnmpEEy 



2gAmpV^P\'i-l5cose) 



+gAmp6^\3 - Pcos9) + -fvmp6^P\l + Pcos9) 

{v3,wm) 



(E.17) 



Similar calculations are made for Mh, ' and we finally obtain (97) 

Appendix F 

Nucleoli Form Factors 

To calculate numerical constants that appear in (82) and (83), we use 
the dipole form factors {^p = 1.793, 



F^'\e) - F^''\e) 



F^\k') - F^-\k') 





-1.913): 




1 


(1 


— k'^/niy) 




1 


(1 


— k'^/rriyf'' 




Hp "H /^n 


(1 


— k^/niy)'^ ' 







{1-k-^/mlf' 

where we adopt my = 0.84 GeV. We also adopt the dipole form factors 



(F.l) 
(F.2) 
(F.3) 
(F.4) 



(1 - fcVmf.) 



2 \2 ' 



FA{k^) 



(1- P/m^)^' 



(F.5) 
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for the vector and axial-vector vertices; we take niA = 1-05 GeV. Weak 
magnetism should also be endowed with the form factor 



The integrals of (82) are inverted as 

1 /• d^fc 1 1 



(F.6) 



167r2 



4-^^ = \ J 



mi 



(27r)4 (A;2)2 {p^ - - ml 

x{F[^\k')+F!f\k')]Fy{k'),{F.l) 



1 f d'^k 1 



1 



(27r)4 (A;2)2 - fc)2 - mj 



{F[^\k^) + F^^\k^)]Fv{k^ 



(F.8) 



We use integral representations employing the Feynman trick to compute 
numerically the constants in (82) and (83), 



J (2^ 



d^k 1 



(27r)4 (fc2)2 (p2 - kf - m2 (1 - A;2/A2)2 (i _ A;2/A2)2 

xle-^{k.pA 



mt. 



AgA?/(Ao,Ai) + ^^^^{Ag/(Ao,0)-A?/(0,Ai)} , (F.9) 



where 



/(Ao,Ai) = dx dyxil-x-y) 

(2-6) 



|y2^2 _|. r^p^2 _|_ ^1 _ 3, _ y)A2|' 

(4^0 - 2)y2m2 



|y2m2 + xAg + (1 - X - y)Af | 
and 6 is either 1/2 or 2. A similar expression is used to evaluate (98). 



(F.IO) 
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The evaluation of (99) goes similarly, starting from 

i cM ^ ± [ d^k {k-p2) 1 
167r2 mjj (27r)4 k'^ (p^ - k)^ - 

xFw{e){FiP\k')+Fi^\k')}. (F.ll) 
Our numerical computation employs the Feynman integral 



(27r)4 fc2 (p2-A:)2-m2(l-F/A2)2(i_A;2/A2)2 



(F.12) 



{/i(l,Ao,Ai)-/i(0,Ao,Ai)}, (F.13) 



where 



^(C,Ao,Ai) = dx dyx{l-x-y) 

x-f — T3- 

[y^^^mj + arA2 + (1 - X - y)A2| 

Likewise, S^"'^ is computed via the same type of parameter integrals. 
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